KERNEL FEATURE CROSS-CORRELATION FOR UNSUPERVISED QUANTIFICATION OF DAMAGE FROM WINDTHROW IN FORESTS

In this study estimation of tree damage from a windthrow event using feature detection on RGB high resolution imagery is assessed. An accurate quantitative assessment of the damage in terms of volume is important and can be done by ground sampling, which is notably expensive and time-consuming, or by manual interpretation and analyses of aerial images. This latter manual method also requires an expert operator investing time to manually detect damaged trees and apply relation functions between measures and volume which are also error-prone. In the proposed method RGB images with 0.2 m ground sample distance are analysed using an adaptive template matching method. Ten images corresponding to ten separate study areas are tested. A 13 7 13 pixels kernel with a simplified lin ear-feature representation of a cylinder is applied at different rotation angles (from 0\ub0 to 170\ub0 at 10\ub0 steps). The higher values of the normalized cross-correlation (NCC) of all angles are recorded for each pixel for each image. Several features are tested: percentiles (75, 80, 85, 90, 95, 99, max) and sum and number of pixels with NCC above 0.55. Three regression methods are tested, multiple regression (mr), support vector machines (SVM) with linear kernel and random forests. The first two methods gave the best results. The ground-truth was acquired by ground sampling, and total volumes of damaged trees are estimated for each of the 10 areas. Damaged volumes in the ten areas range from 3c1.8 7 102 m3 to 3c1.2 7 104 m3. Regression results show that smv regression method over the sum gives an R-squared of 0.92, a mean of absolute errors (MAE) of 255 m3 and a relative absolute error (RAE) of 34% using leave-one-out cross validation from the 10 observations. These initial results are encouraging and support further investigations on more finely tuned kernel template metrics to define an unsupervised image analysis process to automatically assess forest damage from windthrow

wkbj analysis in the periodic wake of a cylinder

Abstract The nature of the three-dimensional transition arising in the flow past a cylinder is investigated by applying the Lifschitz–Hameiri theory along special Lagrangian trajectories existing in its wake. Results show that the von Karman street is unstable with regard to short-wavelength perturbations. The asymptotic analysis predicts the possible existence of both synchronous (as modes A and B) and asynchronous (as mode C) instabilities, each associated to specific Lagrangian orbits. The proposed study provides useful qualitative information on the origin of the different modes but no quantitative agreement between the growth rates predicted by the asymptotic analysis and by a global stability analysis is observed. The reasons for such mismatch are briefly discussed and possible improvements to the present analysis are suggested

Analysis of unmitigated large break loss of coolant accidents using MELCOR code

In the framework of severe accident research activity developed by ENEA, a MELCOR nodalization of a generic Pressurized Water Reactor of 900 MWe has been developed. The aim of this paper is to present the analysis of MELCOR code calculations concerning two independent unmitigated large break loss of coolant accident transients, occurring in the cited type of reactor. In particular, the analysis and comparison between the transients initiated by an unmitigated double-ended cold leg rupture and an unmitigated double-ended hot leg rupture in the loop 1 of the primary cooling system is presented herein. This activity has been performed focusing specifically on the in-vessel phenomenology that characterizes this kind of accidents. The analysis of the thermal-hydraulic transient phenomena and the core degradation phenomena is therefore here presented. The analysis of the calculated data shows the capability of the code to reproduce the phenomena typical of these transients and permits their phenomenological study. A first sequence of main events is here presented and shows that the cold leg break transient results faster than the hot leg break transient because of the position of the break. Further analyses are in progress to quantitatively assess the results of the code nodalization for accident management strategy definition and fission product source term evaluation

Pool temperature stratification analysis in CIRCE-ICE facility with RELAP5-3D© model and comparison with experimental tests

In the frame of heavy liquid metal (HLM) technology development, CIRCE pool facility at ENEA/Brasimone Research Center was updated by installing ICE (Integral Circulation Experiments) test section which simulates the thermal behavior of a primary system in a HLM cooled pool reactor. The experimental campaign led to the characterization of mixed convection and thermal stratification in a HLM pool in safety relevant conditions and to the distribution of experimental data for the validation of CFD and system codes. For this purpose, several thermocouples were installed into the pool using 4 vertical supports in different circumferential position for a total of 119 thermocouples [1][2]. The aim of this work is to investigate the capability of the system code RELAP5-3D (c) to simulate mixed convection and thermal stratification phenomena in a HLM pool in steady state conditions by comparing code results with experimental data. The pool has been simulated by a 3D component divided into 1728 volumes, 119 of which are centered in the exact position of the thermocouples. Three dimensional model of the pool is completed with a mono-dimensional nodalization of the primary main flow path. The results obtained by code simulations are compared with a steady state condition carried out in the experimental campaign. Results of axial, radial and azimuthal temperature profile into the pool are in agreement with the available experimental data Furthermore the code is able to well simulate operating conditions into the main flow path of the test section

Pre-test analysis of protected loss of primary pump transients in CIRCE-HERO facility

In the frame of LEADER project (Lead-cooled European Advanced Demonstration Reactor), a new configuration of the steam generator for ALFRED (Advanced Lead Fast Reactor European Demonstrator) was proposed. The new concept is a super-heated steam generator, double wall bayonet tube type with leakage monitoring [1]. In order to support the new steam generator concept, in the framework of Horizon 2020 SESAME project (thermal hydraulics Simulations and Experiments for the Safety Assessment of MEtal cooled reactors), the ENEA CIRCE pool facility will be refurbished to host the HERO (Heavy liquid mEtal pRessurized water cOoled tubes) test section to investigate a bundle of seven full scale bayonet tubes in ALFRED-like thermal hydraulics conditions. The aim of this work is to verify thermofluid dynamic performance of HERO during the transition from nominal to natural circulation condition. The simulations have been performed with RELAP5-3D© by using the validated geometrical model of the previous CIRCE-ICE test section [2], in which the preceding heat exchanger has been replaced by the new bayonet bundle model. Several calculations have been carried out to identify thermal hydraulics performance in different steady state conditions. The previous calculations represent the starting points of transient tests aimed at investigating the operation in natural circulation. The transient tests consist of the protected loss of primary pump, obtained by reducing feed-water mass flow to simulate the activation of DHR (Decay Heat Removal) system, and of the loss of DHR function in hot conditions, where feed-water mass flow rate is absent. According to simulations, in nominal conditions, HERO bayonet bundle offers excellent thermal hydraulic behavior and, moreover, it allows the operation in natural circulation

Post-test simulation of a PLOFA transient test in the CIRCE-HERO facility

CIRCE is a lead–bismuth eutectic alloy (LBE) pool facility aimed to simulate the primary system of a heavy liquid metal (HLM) cooled pool-type fast reactor. The experimental facility was implemented with a new test section, called HERO (Heavy liquid mEtal pRessurized water cOoled tubes), which consists of a steam generator composed of seven double-wall bayonet tubes (DWBT) with an active length of six meters. The experimental campaign aims to investigate HERO behavior, which is representative of the tubes that will compose ALFRED SG. In the framework of the Horizon 2020 SESAME project, a transient test was selected for the realization of a validation benchmark. The test consists of a protected loss of flow accident (PLOFA) simulating the shutdown of primary pumps, the reactor scram and the activation of the DHR system. A RELAP5-3D© nodalization scheme was developed in the pre-test phase at DIAEE of “Sapienza” University of Rome, providing useful information to the experimentalists. The model consisted to a mono-dimensional scheme of the primary flow path and the SG secondary side, and a multi-dimensional component simulating the large LBE pool. The analysis of experimental data, provided by ENEA, has suggested to improve the thermal–hydraulic model with a more detailed nodalization scheme of the secondary loop, looking to reproduce the asymmetries observed on the DWBTs operation. The paper summarizes the post-test activity performed in the frame of the H2020 SESAME project as a contribution of the benchmark activity, highlighting a global agreement between simulations and experiment for all the primary circuit physical quantities monitored. Then, the attention is focused on the secondary system operation, where uncertainties related to the boundary conditions affect the computational results

Opportunistic rain rate estimation from measurements of satellite downlink attenuation: A survey

Recent years have witnessed a growing interest in techniques and systems for rainfall surveillance on regional scale, with increasingly stringent requirements in terms of the following: (i) accuracy of rainfall rate measurements, (ii) adequate density of sensors over the territory, (iii) space‐time continuity and completeness of data and (iv) capability to elaborate rainfall maps in near real time. The devices deployed to monitor the precipitation fields are traditionally networks of rain gauges distributed throughout the territory, along with weather radars and satellite remote sensors operating in the optical or infrared band, none of which, however, are suitable for full compliance to all of the requirements cited above. More recently, a different approach to rain rate estimation techniques has been proposed and investigated, based on the measurement of the attenuation induced by rain on signals of pre‐existing radio networks either in terrestrial links, e.g., the backhaul connections in cellular networks, or in satellite‐to‐earth links and, among the latter, notably those between geostationary broadcast satellites and domestic subscriber terminals in the Ku and Ka bands. Knowledge of the above rain‐induced attenuation permits the retrieval of the corresponding rain intensity provided that a number of meteorological and geometric parameters are known and ultimately permits estimating the rain rate locally at the receiver site. In this survey paper, we specifically focus on such a type of “opportunistic” systems for rain field monitoring, which appear very promising in view of the wide diffusion over the territory of low‐cost domestic terminals for the reception of satellite signals, prospectively allowing for a considerable geographical capillarity in the distribution of sensors, at least in more densely populated areas. The purpose of the paper is to present a broad albeit synthetic overview of the numerous issues inherent in the above rain monitoring approach, along with a number of solutions and algorithms proposed in the literature in recent years, and ultimately to provide an exhaustive account of the current state of the art. Initially, the main relevant aspects of the satellite link are reviewed, including those related to satellite dynamics, frequency bands, signal formats, propagation channel and radio link geometry, all of which have a role in rainfall rate estimation algorithms. We discuss the impact of all these factors on rain estimation accuracy while also highlighting the substantial differences inherent in this approach in comparison with traditional rain monitoring techniques. We also review the basic formulas relating rain rate intensity to a variation of the received signal level or of the signal‐to-noise ratio. Furthermore, we present a comprehensive literature survey of the main research issues for the aforementioned scenario and provide a brief outline of the algorithms proposed for their solution, highlighting their points of strength and weakness. The paper includes an extensive list of bibliographic references from which the material presented herein was taken

Novel Approach to Rainfall Rate Estimation based on Fusing Measurements from Terrestrial Microwave and Satellite Links

Reliable, cheap and accurate measurements of rainfall rate are growing to be more and more important in many sectors as: meteorology, agriculture, flood warning, and weather forecasting. Recently, indeed, the development of novel competitive techniques has been pushed in order to improve accuracy and reliability performance, such as commercial microwave links and broadcast satellite links. The aim of the current paper is to extend previous works of the literature based on land wireless links only. The basic idea consists in synergically employing both land and satellite based approaches together, by collecting and properly fusing the corresponding measurements. To this end, an iterative optimization procedure has been developed. As shown by numerical results, the proposed procedure gives the estimated rainfall map with a considerable accuracy and improved performance respect to the conventional algorithm based on terrestrial link only

Centre-Manifold Reduction of Bifurcating Flows

In this paper we describe a general and systematic approach to the centre-manifold reduction and normal form computation of flows undergoing complicated bifurcations. The proposed algorithm is based on the theoretical work of Coullet & Spiegel (SIAM J. Appl. Maths, vol. 43(4), 1983, pp. 776821) and can be used to approximate centre manifolds of arbitrary dimension for large-scale dynamical systems depending on a scalar parameter. Compared with the classical multiple-scale technique frequently employed in hydrodynamic stability, the proposed method can be coded in a rather general way without any need to resort to the introduction and tuning of additional time scales. The method is applied to the dynamical system described by the incompressible NavierStokes equations showing that high-order, weakly nonlinear models of bifurcating flows can be derived automatically, even for multiple codimension bifurcations. We first validate the method on the primary Hopf bifurcation of the flow past a circular cylinderand after we illustrate its application to a codimension-two bifurcation arising in the flow past two side-by-side circular cylinder

Boundary-layer Flows Past an Hemispherical Roughness Element: DNS, Global Stability and Sensitivity Analysis

We investigate the full three-dimensional instability mechanism arising in the wake of an hemispherical roughness element immersed in a laminar Blasius boundary layer. The inherent three-dimensional flow pattern beyond the critical Reynolds number is characterized by coherent vortical structures called hairpin vortices. Direct numerical simulation is used to analyze the formation and the shedding of hairpin packets inside the shear layer. The first bifurcation characteristics are investigated by global stability tools. We show the spatial structure of the linear direct and adjoint global eigenmodes of the linearized Navier-Stokes operator and use structural sensitivity analysis to locate the region where the instability mechanism acts. Results show that the “wavemaker” driving the self-sustained instability is located in the region immediately past the roughness element, in the shear layer separating the outer flow from the wake region